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i : 

Wc investigate the dynamo problem in the limit of small magnetic Prandtl number (Pm) using a 
jy^hell model of magnetohydrodynamic turbulence. The model is designed to satisfy conservation laws 
^Jpf total energy, cross helicity and magnetic helicity in the limit of inviscid fluid and null magnetic 
• ^Hdiffusivity. The forcing is chosen to have a constant injection rate of energy and no injection of 
^-kinetic helicity nor cross helicity. We find that the value of the critical magnetic Reynolds number 
^Tjflm) saturates in the limit of small Pm. Above the dynamo threshold we study the saturated regime 
^J^versus Rm and Pm. In the case of equipartition, we find Kolmogorov spectra for both kinetic and 
C^jjiagnetic energy except for wave numbers just below the resistive scale. Finally the ratio of both 
iissipation scales (viscous to resistive) evolves as Pm" 3 ' 4 for Pm < 1. 



Introduction 
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in ; 

^""jfvlost of astrophysical bodies possess or have had in their history their own 
^magnetic fields. In most cases their generation rely on inductive processes 
reproduced by the turbulent motion of the electroconducting fluid within the 
^body [TJ. An important parameter of the problem is the magnetic Prandtl 

^number defined by Pm = vjr] where v is the viscosity and ij the magnetic 
• ^ndiffusivity of the fluid. In the "magnetic" universe Pm varies from values as 

^large as 10 14 for the interstellar medium [2] to values as small as 10~ 6 for 
,i5tlie iron core of planets or stellar plasmas. This large spectrum of possible 

QPm values implies strong differences between possible generation mechanisms. 

t>In some sense Pm is a measure of the kinetic energy spectrum available for 
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generating magnetic energy. When Pm > 1 the resistive scale is smaller than 
the viscous scale implying that all velocity scales are available for generating 
some magnetic field. In the other hand for Pm < 1, only the velocity scales 
larger than the resistive scale are available for the magnetic field generation. In 
that case, the velocity scales smaller than the resistive scale are enslaved to the 
larger scales and in essence they stay passive in the generation process. Besides 
this is why the large eddy simulation technique may be recommended in that 
case [3]. Therefore, at first sight one can expect that dynamo action is all 
the more difficult to obtain since Pm is smaller in reason of a smaller velocity 
spectrum available for the magnetic generation. This is indeed what comes out 
from recent numerical simulations [U El El EE] ( see a ^ so anc ^ re f erences 
therein for an alternative approach). Though, we have evidence of magnetic 
field in planets and stars, and dynamo action has also been reproduced in 
experiments working with liquid sodium for which Pm is small (~ 1(T 6 ) PU 
[T2|I13U14] . These experiments and further devices in preparation [150161117] are 
designed in such a way that the dynamo mechanism is produced by the large 
scale of the flow due to an appropriate large scale forcing. The turbulence 
naturally developing at smaller scales may play a role though this is still 
unclear [181 119|. [3l I20| . I21j . In these experiments, the choice of the forcing is 
based on the hypothesis that it is the stationary part of the large scale flow 
which should be important for the generation mechanism. A number of flow 
geometries studied in the past turned out to be good candidates for such 
experiments [2211231123]. 

In the present paper we are interested in the possibility for a Kolmogorov 
type turbulent flow to generate dynamo action at low Pm, without need for 
a large scale motion controlling the generation mechanism. We expect the 
eddies having the highest shearing rate to be the more active for generating the 
magnetic field, at least during the kinematic stage of magnetic field growth. As 
in Kolmogorov turbulence ui/l ~ Z~ 2//3 , these eddies correspond to the smallest 
available scale which is the viscous scale for Pm > 1 [25] and the resistive scale 
for Pm < 1 [9]. Eventually the magnetic field will then spread out to larger 
scales due to the nonlinear interactions. This problem is hard to solve by 
direct numerical simulation for it needs high resolution in order to describe 
magnetic phenomena adequately [26]. Some results have been obtained using 
the EDQNM closure applied to the MHD equations [27] near the critical Rm 
and for arbitrary low values of Pm. Here we want to investigate arbitrary large 
values of Rm and small values of Pm. For that we use a shell model of MHD 
turbulence introduced by Frick and Sokoloff [28] . This model is the successor 
of several other shell models for MHD turbulence [29,30,31,32,33,34,35] but it 
is the only one to conserve all integrals of motions including magnetic helicity 
(or kinetic helicity for the non magnetic case). It is based on the so-called 
GOY hydrodynamic shell model [36,37,38,39]. In [28], Frick and Sokoloff have 
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derived a model which represents either 2D or 3D MHD turbulence, depending 
on the choice of two parameters. As in real MHD turbulence the 2D model 
leads to the impossibility of dynamo action |40j . This shows that in spite that 
such a shell model is a drastic simplification of the real MHD turbulence, 
ignoring for example the geometrical structures of the motion and magnetic 
field, it contains enough features to make the difference between the 2D and 3D 
problems (see also [Hj). It also reproduces quite well the structure functions at 
different orders of real MHD turbulence. Here we consider only the 3D model 
herein after referred to as FS98. This model has also been used by Lozhkin et 
al. [32] to show that small scale dynamo is possible at low Pm, contrary to 
the hypothesis put forward by Batchelor [33]. 

Giulani and Carbone [H] have shown that long runs with the FS98 model 
lead inevitably towards a "dynamical alignement" stopping the nonlinear 
transfer towards the smaller scales. Giulani and Carbone [H] suggested that 
this problem might be overcome with an other choice of the external driving 
force. This is what we have done here, adopting a forcing in such a way that it 
acts on several scales and depends on time with a random phase at each forc- 
ing scale (see section [231) . Finally, we took care to have long runs well beyond 
any transient state, in order to have good statistics and reliable results. 



2 Shell model for MHD turbulence 



2.1 Model equations 



The shell model is built up by truncation of the Navier-Stokes and induction 
equations. We define logarithmic shells, each shell being characterized by one 
real wave number k n = ko\ n and dynamical complex quantities U n and B n 
representative of the velocity and magnetic fluctuations for wave vectors of 
norm ranging between k n and k n +\. The parameter A is taken equal to the 
gold number (1 + \/5)/2 for it optimizes the resolution [33]. The model is 
described by the following set of equations (0 < n < N) 

d t U n = ik n (Q n (U, U, a) - Q n (B, B, a)) 

- uk 2 n U n + F n , (1) 
d t B n = ik n (Q n (U, B, b) - Q n (B, U, b)) 

- r]k 2 n B n , (2) 

where 

Q n (X, Y, c) = Cl X* n+1 Y* +2 + c 2 x:_ x y n * +1 
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+ c 3 X*_ 2 Y*_ v (3) 

represents the nonlinear transfer rates with the four neighbouring shells n — 2, 
n — 1, n + 1 and n + 2. In addition we have to take U-2 = U-\ = Un+i = 
Un+2 = and _B_2 = B-i = -B/v+i = B^+2 = 0. The parameter F n is 
the forcing at shell n. The time unit is defined by the turnover time of the 
largest scale r = (|f7o|fco) _1 - To determine the complex coefficients dj and bj, 
j = 1, 2, 3 we apply the property that the total energy E to t, cross-helicity Tic 
and magnetic helicity TCb must be conserved in the limit of non- viscous and 
non-resistive limit v = rj = 0. In our shell model, these quadratic quantities 
write in the following form 

TV 

I 

Etot 



^(\U n \ 2 + \B n f), (4) 

n=0 
1 - 

n c = -^2(U n B* n + B n U:), (5) 

n=0 
1 N 



2 

n=0 



leading to a\ = 1, a 2 = (1 - A)A~ 2 , a 3 = -A -3 , h = b 2 = b 3 = (A(l + A)) -1 . 
In the pure hydrodynamic case (B n = 0) the original GOY model is recovered 
satisfying, in addition to ([!]), the conservation of the kinetic helicity 



1 - 

Hu = -Y,(-mUn\ 2 k n . (7) 



n=0 



2.2 Forcing and initial conditions 

The forcing is chosen in order to control the injection rate of kinetic energy, 
cross and kinetic helicities. For that we spread the forcing on three neighbour- 
ing shells nf,Tif + 1 and n/ + 2 with F nf +j = fje l< ^ j , j = 0, 1, 2 where the fj 
are positive real quantities and where the (f>j S [0, 2tt] are random phases. In 
that case the forcing is ^-correlated. Alternatively we also used a forcing for 
which the phases 4>j are constant during a certain time r c , which can be inter- 
preted as a finite correlation time. In fact this does not make much difference 
either on the autocorrelation functions of U n nor on the subsequent results. 
Therefore it is sufficient to use random phases. As we are interested to inject 
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neither kinetic helicity nor cross-helicity, the forcing functions must satisfy 

n f +2 

- £ KFn + U n F* = e, (8) 

n=n } 

n f +2 

J2(-l) n k n (U*F n + U n F*) = 0, (9) 

n=nf 

n f +2 

B*F n + B n F* = 0. (10) 

n=n } 

where e is the rate of kinetic energy supplied to the system. Therefore for a 
given set of random cfrj (j = 0, 1, 2), the fj depend on the Uj and Bj (j = 0, 1, 2) 
which expressions are given in Appendix. For some arbitrary initial conditions 
on Uj (j = 0, 1) of small intensity (~ 10 -6 ) we let the hydrodynamic evolve 
until it reaches some statistically stationary state. Then introducing at a given 
time some arbitrary non zero values of Bj (j = 0, 1) of small intensity (~ 1CP 6 ) 
we solve the full problem until a statistically stationary MHD state is reached. 
The time of integration needed to obtain good statistics depends on v and r\ 
but typically it is equal to several hundreds of the large scale turn-over time. 



2.3 Input and output 

The input parameters of the problem are v, rj, the forcing shell rif, the rate of 
injected kinetic energy e and the number of shells N. In the rest of the paper 
we take e = 1. 

As output we define the kinetic and magnetic energy for the shell n by 

E U (n)= l -\U n \ 2 and E B (n) = l -\B n \ 2 , (11) 
the total kinetic and total magnetic energy by 

N N 

Eu = ^2E u (n) and E B = ^E B (n) (12) 
n=0 n=0 

and the total energy by 



Etot = E\j + Eb ■ 



(13) 
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Following we define the spectral energy fluxes from the inside of the [/(or 
-B)-sphere (shells with k < k n ) to the outside of the [/(or l?)-sphere (shells 
with k > k n ). We note for example 11^5 (n) the energy flux from the inside of 
the B-sphere to the outside of the [/-sphere. Then we have 

71- 1 

Y,^{k 3 U*Q 3 {U,U,a)} (14) 

3=0 

71-1 

^{-kjUjQjiB, B, a)} (15) 

3=0 

71-1 

Y Qi-kjB^QjiB, U, b)} (16) 

3=0 

71-1 

^{kjB^Q^B,^}. (17) 

3=0 

In FS98 the time average of II^(n) is denoted U n . We also define the energy 
fluxes from the inside of the [/-and--B-spheres to the outside of the [/-sphere 
or B-sphere by 

n u (n)=nV<(n) + Il§<(n) (18) 
n B (n)=n^<(n)+n|<(n) (19) 

and the total energy flux by 



II^<(n) = 
U§<(n) = 
U u B <(n) = 
Hf<(n) = 



U tot (n)=U u {n)+U B (n). (20) 

We define the viscous and resistive dissipation rates D u (n) and D B (n) in shell 

n, by 

D u (n) = vk 2 n \U n \ 2 (21) 
D B {n) = ^kl\B n \ 2 (22) 



and the total dissipation rate by 



N 

D tot = Y J {D U {n) + D B {n)). 

n=0 



(23) 
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With these definitions we obtain the following shell-by-shell energy budget 
equations: 

n n 

dt E u (j) + Uu(n) = -J2 D u (j) + e (24) 

j=0 j=0 
n n 

d t Y J E B { 3 )+Ii B {n) = -Y J D B {j). (25) 

j=0 j=0 

For a statistical stationary solution (dt(E u (j)) = dt(E B (j)) = 0) we have 
then 

n n 

(Vtot(n)) = - *£(D u (j)) - Y,(D B U)) + e. (26) 

3=0 j=0 

where here and after ( ) denotes time averaged quantities. 
We define the kinetic and magnetic Reynolds numbers as 

Re = {E tot ) 2 /(u(D tot )) (27) 
Rm=(E w ) 2 /(# w ». (28) 

Finally, following [47], we define the viscous (resp. resistive) scale k~ l (resp. 
k~ l ) as the one at which the viscous (resp. Ohmic) decay time t v = (vk^)^ 1 
(resp. r n = (rjkn)" 1 ) becomes comparable to the typical turn-over time tjj = 

(KiPnl 2 ) 1 / 2 )- 1 . 



3 Hydrodynamics 

Choosing the appropriate forcing corresponding to B n = we present in FigQ] 
some results concerning the pure hydrodynamic case for v = 10 -8 and nj = 8. 
In this case the forcing is <5-correlated. Though the autocorrelation function, 
defined by 

( , _ JU*(t)Un(t + T) + U n (t)U*(t + T)dt 

2^ J U*{t)U n (t)dt f U n {t + T)U*(t + T)dt 

and plotted in FigfjJ,, is far from being the one of a <5-correlated velocity 
contrary to the Kasantzev model [9j. We also made comparisons with a finite 
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Figure 1. Hydrodynamic case for v = 10 -8 and a forcing scale (arrow) corresponding to nj = 8. 
The output Reynolds number is Re = 8 10 7 . In (a), the autocorrelation function cor(n,r) for a 
5-correlated forcing is plotted versus r and for several shells n. In (c), the turn-over (black dots) 
and dissipation (straight line) characteristic times are plotted versus log 10 k. In (b), the energy 
spectrum is plotted versus log 10 k and the A; -2 / 3 slope (full line) is plotted for comparison. In (d), 
the energy flux (black dots) and the dissipation 5~Jj=o D (j) (gray dots) are plotted versus log 10 k. 



correlation time forcing without finding any significant differences. Therefore 
the (^-correlated forcing does not seem to be an issue in our problem. 

The kinetic energy spectrum (Fig(Tb, black dots) of the stationary statistical 
state is found to be in /c -2 / 3 (which corresponds to a Fourier energy spectrum 
of fc _5//3 as expected in Kolmogorov turbulence). In FigQJ, the spectral flux 
Hjj{n) (black dots) and the dissipation ^™ =0 D u (j) (gray dots) are found to 
satisfy the kinetic energy budget (|24p with e = 1. In addition, in the inertial 
range we find that Hu(n) ~ e and Yl^o^U) ~ as predicted by a Kol- 
mogorov turbulence. After the viscous scale, II;y(n) ~ and /Cj=o D (j) ~ e. 

As previously defined, the viscous scale is the one at which the viscous 
decay time r v = (z^/c 2 ) -1 (full curve of FigfTb) becomes comparable to the 
typical turn-over time ly = {k n U n )~ l (black dots of FigfTb). This leads to 
k u ~ 10 6 and compares indeed very well with the Kolmogorov dissipation 
scale k^ 1 ~ (i^/e) 1 / 4 . Finally the little bump of Uu(n) (black dots FigfTbl) 
just before the viscous scale looks like a bottle-neck effect [38]. 
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Figure 2. Quadratic quantities (a) E to t, (b) E v , (c) E B , (d) Hu, (e) H B and (f) Hc/V^uEb 
versus time, for v = 10~ 9 and n = 10 -6 . 



4 Dynamo action 



4.1 Time evolution of quadratic quantities 

Here we start with a typical example of magnetic generation for v = 10 -9 and 
7] = 1CP 6 (Pm = 10~ 3 ). In Fig. [2] the different quadratic quantities defined 
in dU), ([5|), ©, ([7|) and (fT2|) are plotted versus time. A coarse time sampling 
has been chosen here for a better representation of the results and is not rele- 
vant of the actual time step used for the numerical calculations. The kinetic, 
magnetic and total energies have reached a statistical stationary steady state 
after a few hundred time steps. The fluctuations of these quantities are quite 
important due to the small values of v and rj. The kinetic helicity though 
its average is close to zero, shows strong fluctuations. In the other hand the 
magnetic helicity stays very small. Finally the relative cross helicity defined 
by Ti-c/ V^i/Eb oscillates around zero. The fact that this latter quantity does 
not reach an asymptotic limit of ±1 shows that there is no "dynamical aligne- 
ment". Therefore we are confident that our choice of forcing overcomes the 
problem raised by Giulani and Carbone |41j . 



4.2 Spectrum analysis 

In Fig. [3] we show the kinetic and magnetic spectrum at four successive times 
for again v = 10~ 9 and rj = 10~ 6 (Pm = 10~ 3 ). Each snapshot corresponds 
to an average over a not so large amount of time which explains why at early 
time the kinetic spectrum is not very smooth at large scales. In the early time, 
when the magnetic field is still not significant, the kinetic energy spectrum has 
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Figure 3. Kinetic (black dots) and magnetic (gray dots) spectra at four successive times (from (a) 
to (d)) for nj = 4, v = 10 -9 and rj = 10 -6 . See also the movie energyl.mpg in which log 10 E u (n) 
and log 10 E B (n) are plotted versus log 10 k with respectively red and blue dots. 



a slope in £;~ 2 / 3 (corresponding to a Fourier spectrum in £;~ 5 / 3 ). Then, as Rm 
is much larger than the critical value of the dynamo instability, the magnetic 
energy starts to grow (Fig. [3h). We expect magnetic energy to be initially 
amplified by the eddies having the highest sharing rate, i.e. the smallest scale 
eddies. As Pm < 1, the smallest eddies available for dynamo action correspond 
to eddies at resistive scale. This is indeed what we find, as here, the resistive 
scale (defined as in section I2T31) corresponds to log 10 k v ~ 4.1. We note that 
the Kolmogorov resistive scale given by k^ ~ (e/77 3 ) 1 / 4 (see section l4~3l) with 
rj = 10~ 6 , leads to a slightly higher value log 10 k v ~ 4.5. 

As Rm is sufficiently large, at subsequent times the magnetic energy reaches 
the level of kinetic energy (Fig. [3b). At that time the kinetic spectrum is 
not influenced yet by the nonlinear feedback of the magnetic field and is still 
in A; -2 / 3 . Then the dynamical equilibrium between the magnetic and velocity 
fields settles down (Fig. [3H) . A striking feature of this equilibrium is the change 
of slope (from -2/3 to ~ -1) of the kinetic energy spectrum for k < k v while 
the magnetic spectrum is slightly above the kinetic spectrum. We also note 
that the viscous dissipation scale has increased (the right part of the kinetic 
spectrum drifting to the left). This probably comes from the fact that there 
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logiok logiok logi k 

Figure 4. Kinetic (black dots) and magnetic (gray dots) spectra for u = 10~ 9 and for Pm = 
(a) 1CJ- 2 , (b) 1CT 1 , (c) 10° and Re = (a) 6.5 10 9 , (b) 4.4 10 9 , (c) 4.4 10 9 . The forcing scale 

corresponds to nj = 4. 




is less energy to dissipate by viscosity than at earlier time because of the 
additional Joule dissipation. 

When changing the value of Pm while keeping the same value of v and cal- 
culating again the final statistically stationary state, we observe again (Fig. [4]) 
a deviation of the kinetic energy slope from -2/3 to ~ -1 whatever the value 
of Pm. To understand better these spectra, we plotted several fluxes in Fig. El 
for v = 1CT 9 and Pm = 1CT 3 . 

Looking at curve (a) which represents the total flux H tot (n) versus log 10 k, 
one can distinguish three plateaus: the first one corresponds to scales larger 
than the resistive scale (1 < log 10 k < 3), the second one for scales smaller 
than the resistive scale but larger than the viscous scale (log 10 k ~ 5), and 
the third one for scales smaller than the viscous scale (log 10 k > 7). The drop 
from the first to the second plateau corresponds to the ohmic dissipation rate 
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e n = Ylj=o D {])■ The drop from the second to the third plateau corresponds 
to the viscous dissipation rate e v = Ylj=o-D U (j)- We clearly have e = e u + e v 
as expected from (j26j) for n = N . 

The curve (b) corresponds to Hjj(n) versus log 10 k with two plateaus, de- 
pending if the scale is larger or smaller than the viscous scale. The first 
plateau (A; < 6) corresponds to IIj/(n) ~ e and the second one {k > 7) to 
II{/(n) ~ e — e v = e v . In particular, there is no clear change of Il(/(n) just 
before the resistive scale that could explain the change of slope of the kinetic 
energy spectrum as previously pointed out. 

Now let us have a look at curve (c). The transfer rate 11^5 (n) is responsible 
for the direct cascade of kinetic energy and would be constant leading to a 
Kolmogorov spectrum if the magnetic field was null (see Fig. [TJ). This would 
remain true for a non zero magnetic field only if the curve (c) was staying flat 
with Il^(n) = e v for 2 < log 10 k < 5.5. In that case the curve (d) would be 
flat as well with 11^ (n) = e v for k > 2. Instead, there is a drop of (n) 
compensated by a symmetric bump of (n) for 2 < log 10 k < 4.5. This drop 
of n^(n) is consistent with a spectrum steeper than A;~ 2 / 3 . Indeed, the bump 
of n^(n) corresponds to some extra energy taken from e and dissipated by 
Joule effect. Then there is less energy to be transferred through the kinetic 
energy cascade. The physical reason why this scenario happens for scales just 
larger than the resistive scale, however is still unclear. 

For the parameters of Fig. [5] the Kolmogorov dissipation scales are given by 
k v = (e/r/ 3 ) 1 / 4 = 10 4 ' 5 and k T] = (e/V 3 ) 1 / 4 = 10 6 ' 75 which correspond quan- 
titatively well with the beginning of the second and third plateau of Tltot (n). 
This shows that the arguments leading to the Kolmogorov dissipation scales 
(see next section) are not affected by the change of spectra slopes observed in 
Fig.H 

Finally for completeness, we produced three movies showing the time evo- 
lution of the spectra of the other quadratic quantities. In u-helicity.mpg, b- 
helicitympg and cross-helicitympg, log 10 Hu(n), log w TtB(n) and log 1Q Hc(n) 
are plotted versus log 10 k where the blue and red dots denote positive and neg- 
ative signs. 



4.3 Dissipation scales ratio 

At the end of section [231 we have already explained how we identify the viscous 
and resistive scales k u and k v , by comparing the turn over time to the respec- 
tive dissipative times. In Fig. [6] we plot the ratio k^/k^ versus Pm < 1 for 
different values of Re. We find that k u /k v ~ Pm _3//4 . To understand why, it is 
sufficient to say that between k v and k u the kinetic energy obeys a Kolmogorov 
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Figure 6. Ratio ky/k^ versus Pm for different values of indicated in the legend. The straight 
line fc -3 / 4 is plotted (dashed line) for comparison. 



spectrum U(k) = e 1 / 3 k 1 / 3 (see Fig. |5J), leading to r^ 1 = kll(k) = e 1 / 3 ^ 3 . 
Comparing r^ 1 with respectively r" 1 = vk 2 and r^ 1 = r\k 2 leads [17] to the 
dissipation scales k v ~ (z^/e) -1 / 4 and fc^ ~ (?7 3 /e) — 1//4 . This in turn leads to 
a dissipation scales ratio in Pm -3 / 4 . 



4.4 Route to saturation 

In this section we study the influence of Pm on the way the dynamo saturates. 
For that we calculate the ratio of magnetic to kinetic energy Eb/Ejj, Eb 
and Ejj being defined as in (fT2j) . In Fig. [7J Eb/Ejj is plotted versus Rm 
for three values of Pm. We note that for Rm much larger than the critical 
value, the level of saturation Eb/Eu may go beyond 1 for Rm ~ 10 5 . Such 
a super saturation state could be expected from the spectra of Fig. [U At 
the threshold, the slope of Eb/Ejj versus Rm follows a turbulent scaling of 
the form Eb/Eu ~ (Rm — Rm c )/Rm^ as expected by Petrelis and Fauve [55] . 
Indeed as in this case the threshold Rm c does not vary very much with Pm, the 
slopes at Rm = Rm c are similar. This is to contrast with the laminar scaling 
Eb/Eu ~ Pm(Rm — Rm c )/Rm^ [49] which would lead to a quasi-horizontal 
slope for Pm = 10 -4 . 
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Figure 7. The energy ratio Eb/Eu versus Rm for rif = 4 and three values of Pm. 

4.5 Dynamo threshold 

In Fig. [8] the dynamo threshold Rm c is plotted versus Pm -1 for nj = 4. For 
increasing values of Pm -1 up to 10 3 the threshold first increases in accordance 
with previous direct numerical simulations [U[5l[6l[71[8l[9] . However, for values 
of Pm -1 larger than 10 3 the threshold Rm c is found to reach a plateau. 

For each value of Pm, the vertical bar around Rm c corresponds to values 
of Rm for which the magnetic solution is erratic. In other words, below the 
bars there is no dynamo action and above the bars there is a well define 
statistically stationary magnetic solution. In between though we do not observe 
intermittency as in |50U51j . the dynamo is irregular, the mean magnetic energy 
increasing and decreasing versus time. 

4.6 Influence of a forcing scale smaller than the resistive dissipation 
scale 

In Fig. [9l the kinetic and magnetic spectra are plotted for a forcing scale 
smaller than the resistive scale k v . In that case the inertial range does not play 
a role in the magnetic generation and a kinetic spectra in £;~ 2 / 3 is recovered. 

5 Discussion 

In this paper we investigated the fully developed MHD turbulence at magnetic 
Prandtl number lower than unity, using a shell model of MHD turbulence 
with an appropriate forcing. The main results are: 
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Figure 8. Dynamo threshold Rm c versus Pm 1 for nf = 4. 
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Figure 9. Kinetic (black dots) and magnetic (gray dots) stationary spectra for u = 10~ 9 , 
Pm = 10 -7 and a forcing scale corresponding to n/ = 12. See also the movie encrgy2.mpg in which 
log 10 E u (n) and log 10 E B (n) are plotted versus log 10 k with respectively red and blue dots. 



1. For strong MHD turbulent dynamo states (large Rm) we find kinetic and 
magnetic energy spectra close to the Kolmogorov spectrum A; -2 / 3 except 
at scales just larger than the resistive dissipation scale for which there is a 
weaker (stronger) slope of the kinetic (magnetic) spectrum. This corresponds 
to the work of the Lorentz forces which increases with k up to k = k v . 

2. The evaluation of the viscous and resistive dissipation scales are consistent 
with Kolmogorov estimates leading to k u /k n ~ Pm~ 3 / 4 . 

3. At the dynamo threshold Rm c , the ratio of magnetic to kinetic energy scales 
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like Eb/Ejj ~ (Rm — Rm c )/Rm 2 , as predicted by a turbulent scaling [39]. 
4. At very low values of Pm, the dynamo threshold Rm c reaches a plateau. 

Of course all these results rely on the assumption that the interactions 
between the different scales of motion and magnetic field are local interactions, 
each shell interacting with a few shells above and below. We believe that this 
should not make much difference as long as Pm is small, the Kolmogorov 
turbulence being governed by local interaction. In the other hand our results 
can not being tested against the Iroshnikov-Kraichnan A; -3 / 2 Fourier spectrum 
prediction [52] resulting from non local interactions between the flow and 
some large scale magnetic field which could result for example from dynamo 
action. By the way we believe that the fc -3 / 2 slope in FS98 is due to a lack of 
statistics as can be seen from the energy fluxes which are not flat and from 
the corresponding small range of scales. Adding some non local interaction 
with a large scale magnetic field in a local shell model, Biskamp [34] found a 
£,-3/2 si p e) though taking only one such a non local interaction is somewhat 
artificial. Recently Verma [53] revisited the Iroshnikov-Kraichnan theory in 
which he shows that the large scale magnetic field becomes renormalized 
due to the nonlinear term, leading back to the Kolmogorov spectrum. This 
emphasizes the need for a complete nonlocal shell model in which any shell 
could interact with the others. This could be a good test against one theory 
or the other. Such a model would be also welcome for simulations at large 
Pm. Indeed at large Pm we expect the more energetic scales of the flow, 
corresponding to scales close to the viscous scales, to interact directly with the 
smaller scales of the magnetic field. Our local shell model can not catch such 
features and this is why we did not show results at large Pm for they surely 
lack physical ground. A further issue that could be addressed by a nonlocal 
shell model could be to distinguish between a large scale field generated by 
a small scale velocity field resulting from non local interactions (developed 
in the mean field formalism) and a large scale field generated by an "inverse 
cascade" as for example in Fig. [3] or in [53], resulting from local interactions. 
Concerning our local model, we believe that the results presented in Fig. [8] 
showing that the dynamo threshold does not depend on Pm at low values of 
Pm would stay qualitatively the same if additional nonlocal interactions were 
included in the model. Indeed the dynamo threshold corresponds to the growth 
start of the magnetic field which is then still not significant. Therefore any non 
local interactions (e.g. Alfven sweeping effect) might not change the threshold. 
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7 Appendix 

For the pure hydrodynamic case (B n = 0), only the two first conditions ([S]) 
and © are necessary to derive the forcing equations. In that case the forcing 
set writes 

f° = rYXn ^ u T ( 30 ) 

(A + l)u cos(0 o -uj ) 

^ (A + l)u\ cos(0i — ui) ^ ^ 

h = 0, (32) 

while for the full MHD case the forcing set is derived from the three conditions 

da), a and ma 



©, © and mj 

— (1 + A)/ = \b 2 ui cos(6 2 - 4> 2 ) cos(0i - oj\) 
+ \ 2 b\u 2 cos(#i - 4>i) cos(0 2 - Uz) (33) 

— (1 + A)/i = b 2 U C0S(6» 2 - 2 ) cos(0 o - u ) 

- \ 2 b u 2 cos(6 l o - 4>o) cos(0 2 - uj 2 ) (34) 
— (1 + A)/ 2 = -hu cos(6»i - 4>i) cos(0 o - w ) 

— Afco^i cos(0q — 0o ) cos(0i — u\) (35) 



where 



j4 = 62^0^1 £08(6*2 — 02) cos(0o — u>o) cos(0i — u>i) 

+ (A - l)6i U U 2 008(6*1 - 0i) COS(0 O - Wo) COS(0 2 - U) 2 ) 

— XboU\u 2 cos(6 ) o — 0o ) cos(0i — u>i) cos(02 — u> 2 ) (36) 
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and where Uj and ujj (resp. bj and 9j) are the complex modulus and argument 
of U nf+j (resp. B nf+j ). 
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